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Abstract. We present XMM-Newton observations of hard X-ray emission from the held of IC443, a supernova 
remnant interacting with a molecular cloud. Fhe hard emission from the held is dominated by 12 isolated sources 
having 2-10 keV flux <; 10 -14 erg cm -2 s _1 . Only a fraction of the sources are expected to be extragalactic or 
stars on statistical grounds, while the others may be associated with the remnant. We have analyzed near-infrared 
K band and also DSS optical data for all of the detected sources, Hnding that six X-ray sources are located in 
a relatively small 15' x 15' region where there is strong 2.2 fim infrared emission, indicating interaction with a 
molecular cloud. The source 1SAX J0618. 0+2227, the brightest in this region (excluding the plerion), is resolved 
with XMM-Newton into two sources, one of which is extended and has a hard power law spectrum photon index 
< 1.5) and shows some indications of spectral line signatures (e.g. Si), while the other is point-like and has a 
featureless spectrum of steeper photon index ~ 2.2. Possible interpretations of some of the discrete sources in 
terms of interaction between the SNR and the molecular cloud are discussed. 

Key words. Acceleration of particles; Radiation mechanisms: non-thermal; ISM: clouds; ISM: cosmic rays; ISM: 
individual object: IC443; ISM: supernova remnants 



1. Introduction 

The middle-aged SNR IC443 (G189. 1+3.0) appears in 
multiwavelength images as an extended (~ 45 arcmin size) 
source of both nonthermal and thermal continuum emis- 
sions with rich line spectra from radio to X-rays. The ex- 
tended component of the thermal soft X-ray emission from 
IC 443 has been observed wit h the Einstein sat ellite and 
the HEAP 1 A2 experiment l|Petre et alJll98Sft. ROSAT 
llAsaoka fc Aschenbacblll994li and ASCA l|Keohane et alJ 
ll997tlKawasaki et al.ll2002j) . On the other hand, the hard 
X- ray emission f rom IC 443 in the 2-10 keV band, detected 
by IWang et"al ] Jl992ii using data of the Ginqa satellite, 
has b een studied in detail only recently. iKeohane et alJ 
using the ASCA X-ray observatory, mapped for 
the first time the hard X-ray emission of IC443, con- 
cluding that most of the 2-10 keV photons came from 
an isolated emitting feature and from the southeast elon- 
gated ridge of hard emission. They also suggested that 
synchrotron emission of > 10 TeV regime electrons 
provides the explanation for these observed features. 
BeppoSAX observations in the hard X-ray band resolved 
the two ASCA sources into two compact sources of hard 
emission up to 10 keV (1SAX J0617. 1+2221 and 1SAX 
J0618. 0+2227) and a harder component up to 100 keV, 
indicating that a synchroton origin is unlikely, and sug- 
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gesting el ectron acceleration by a s low shock in a molecu- 
lar cloud l|Bocchino fc Bvkovl200o|b While recent observa- 
tions of 1SAX J0617.1+22 21 with Chandra llOlbert et alJ 
l200l|) and XMM-Newton (|Bocchino fc Bvkovll200l|) have 
established the plerionic nature of this source, the debate 
is still open about the nature of 1SAX J0618. 0+2227. 

The morphology of IC443 at high X-ray energies is 
therefore expected to be very interesting. It is impor- 
tant to note that IC443 i s associated with the EGRET 
source 3EG J0617+2238 l|Esposito et alJll996|) . but the 
relation between the remnant and the 7-ray source is not 
yet clear. In fact, the 95% error circle of the most recent- 
derived sky position of 3EG J0617+2238, reported by 
lHartman etafl £999), does not include the plerion (see 
Fig. nj, and the other compact X-ray sources, even if the 
extrapolation of the X-ray flux of 1SAX J0617.1+2221 
and 1SAX J0618.0+2227 to the EGR ET regime indicates 
that there could be a direct relation ^Bocchino fc Bvkovl 
1200(1) . High-resolution hard- X-ray studies are therefore of 
great importance to understand this comp licated object, 
also in the light of the recent discovery bv lUchivama et all 
( 2002) of very similar features in the hard X-ray emission 
of the SNR 7 Cygni observed by the ASCA X-ray satellite. 
Compact sources with hard emission spectra in both IC 
443 and in 7 Cygni could be of similar nature e.g. shocked 
molecular clumps as was predicted by the model of Bykov 
et al. (2000). 
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Another interesting class of hard X-ray sources related 
to the SNR could be fast moving fragments of metal-rich 
ejecta interacting with the ambient medium. The isolated 
X-ray emitting ejecta fragments were discover ed in the 
middle-aged Vela SNR llAschenbach et al1ll995h and also 
in type la Tycho SNR fe.g. iDecourchelle et alJl200l|) . A 
model of X-ray line emission from fast moving ejecta frag- 
ments predicts that the fragments interacting with molec- 
ular clouds are the bright est and they could be detected 
from a few kpc distances l)Bvkovll2002h . 

In this paper, we are not disc ussing the IC 4 43 ple rion, 
which was already descri bed by lOlbert et alJ (j200lT) and 
iBocchino fc Bvkovl (|200l|) . while we will address in detail 
1SAX J0618. 0+2227, a source closely correlated with th e 
molecular line emission region teocchino fc Bvkovll2000l) . 
Moreover, the combination of high sensitivity of the PN 
and MOS cameras aboard XMM-Newton to 2-10 keV pho- 
tons with a few arcseconds angular resolution has allowed 
us to study the hard X-ray emission from the whole ex- 
tended supernova remnant IC443 in great detail, thus dis- 
covering an additional extended component and a new set 
of weaker compact sources. Both were undetected in pre- 
vious BeppoSAX and ASCA observations. In the present 
paper using the Cal/PV dedicated pointings on IC443 we 
present the counting and analysis of all the hard sources of 
X-ray emission within the boundary of the IC443 remnant 
with 2-10 keV flux > 5 x 10~ 14 erg cm -2 s^ 1 , while a 
study of the diffuse hard X-ray emission of this remnant 
will be discussed elsewhere. 

2. Observations and Data analysis 

IC443 has been observed a s part of the Cal/P V phase of 
the XMM-Newton satellite ijjansen et al.l200lj) . In this pa- 
per, we have used th e observations liste d in Tabled Data 
fro m the two MPS former et all 1200 ll) cameras and the 
PN intruder et alJ200ll) camera have been used. MOS and 
PN cameras are CCD arrays which collect X-ray photons 
between 0.1 and 15 keV and have a field of view of 30'. The 
pixel size is 1.1" and 4.1" for MOS and PN respectively, 
while the mirror PSF width is 6" - 15" FWHM-HEW. 
The data have been acquired with the medium filter and 
in full image mode, and therefore the temporal resolution 
is low, 2.5 s and 73 ms for MOS and PN, respectively. The 
worse spatial resolution of the PN is compensated by its 
greater sensitivity, on the average 20-30% more than the 
combined two MOS. 

The Standard Analysis System (SAS) software we have 
used (version 5.3) takes care of most of the required event 
screening. However, we have further screened the data to 
eliminate some residual hot pixels and occasional back- 
ground enhancement due to intense flux of soft protons in 
the magnetosphere. In particular, we have extracted the 
background light-curve at energy > 10 keV and we have 
identified time intervals of unusually high count rates (typ- 
ically more than 25 cnt s _1 ) and removed them from sub- 
sequent analysis. Moreover we have also selected events 
with values of the PATTERN keywords between and 12 



Table 1. Cal/PV Observations of IC443 



Obs. 


Coord (J2000) 


J- lis / J- s 


Date 






ksec 




Ptgl 


6 h 17 m 24.2 s +22 d 43 m 13 s 


23/10 


26 Sep 2000 


ptg2 6 


6 h 16 m 11.2 s +22 d 43 m 29 s 


30/21 


25 Sep 2000 


ptg3 6 


6 h 17 m 24.3 s +22 d 26 m 43 s 


32/24 


27 Sep 2000 


ptg4 


6 h 16 m 11.2 s +22 d 19 m 29 s 


30/23 


28 Sep 2000 



a Unscreened and screened exposure times 
6 Two pointings combined 



for MOS and between and 4 for PN. The exposure time 
of the screened observations is also reported in Table ^ 
In the case of ptg2 and ptg3, which are splitted into two 
subpointings each, we have merged the event files of the 
subpointings before continuing with the analysis. 

2.1. Image mosaic 

For each of the pointings listed in Table we have ex- 
tracted PN images with 8" pixel size and exposure images 
at the same resolution. We have defined two energy bands 
in which to extract images, namely 0.5-2 keV (hereafter 
soft band) and 3-10 keV (hard band). The images in each 
band (and the exposure maps) have been aligned using the 
reference pointings directions. The resultant count image 
has been divided by the exposure maps mosaic to cor- 
rect for exposure and vignetting effect. Fig. ^ shows the 
0.5-2 keV mosaic of the PN pointings. Most of the emis- 
sion in this band is of therm al origin, as established by 
lAsaoka fe Aschenbachl l)l994|) based on ROSAT data. 

Fig.glreports the 3-10 keV PN mosaics of IC443. Most 
of the thermal emission associated to IC443 is not present 
in this band, and therefore this band is very well suited for 
the study of hard compact X-ray sources in the vicinity of 
this SNR. 

2.2. Compact sources 
2.2.1. Source detection 

We have used the ewavelet SAS task (v2.4) to de- 
tect sources in the PN images. This task uses the mex- 
ican hat wavelet algorithm to detect both point and ex- 
tended sources provided that the source extent is not too 
large (< l/8th the image size). A review of detection al- 
gorithms for XMM-Newton , inclu ding ewavelet is pre- 
sented bv lValtchanov et al.l (|200l|> . while further informa- 
tion about this ta s k can be found in the SAS User's Guide. 
IValtchanov et alJ l|200lj) have tested ewavelet with sim- 
ulations of both point-like and extended sources, find- 
ing that the algorithm performs very well for point-like 
sources, whereas it tends to miss more extended sources 
than other algorithms (e.g. emldetect). Since we are 
mostly interested in compact hard X-ray sources, we have 
discarded sources detected at large wavelet scales (> 6 
pixel), and therefore this limitation does not affect our 
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Fig. 1. PN mosaic in the soft (0.5-2 keV) energy band, with surface brightness contours at 5 x 10 -4 , 0.001, 0.002, 
and 0.004 cnt s _1 pix -1 (8" x 8" pixel). We have overlaid in bold style the 12 sources detected, the 95% error circle 
of the position of the EGRET sou rce 3EG J0617+2238, and a contour from the H2 1-0 5(1) line emission map at 
2.122 /zm fro mlBurton et alJ (|l990|) . the emission line which is produced by shocked molecular hydrogen according to 
iBurton et all l|l988)) . 



results. A crucial p arameter for source dete ction is the 
detection threshold; IValtchanov et al I ll200lh show that 
ewavelet finds a number of false detections of 10% of 
the number of input sources in an observation of 10 ks 
with cosmic background, and using a threshold of 4<r. To 
further reduce the number of spurious detections, we have 
adopted a threshold of 5a. 

We have run the algorithm in the soft (0.5-2 keV) and 
hard (3-10 keV) energy band, using wavelet scales in the 
range 2-8 pixels. We have detected 88 sources in the soft 
band, and 12 in the 3-10 keV band. We therefore expect 
less than one spurious source in the hard band. The soft- 



band list is dominated by sources representing fluctuations 
in the thermal soft X-ray emission of this remnant. For 
this reason, we have reported in Table |21 only the 3-10 
keV band detections along with their soft counterpart if 
any. Figs. and [21 also show the positions of the detected 
sources. 

Table reports the source position, the wavelet cor- 
relation coefficients which are a measure of source detec- 
tion significance, the wavelet scale of detection and the 
extension parameter, which is equivalent to a FWHM in 
the approximation of a Gaussian Point Spread Function. 
Simulations with point sources only have shown that this 
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Fig. 2. PN Mosaic in t he 3-10 keV energy ban d. Same overlays a s for F ig. □ Src 5 is 1SAX J0617. 1+2221, the 
plerion nebul a, studied in|Olbert et al. 
J0618+2227 ijBocchino & Bvkovll200rl 



20011) and lBocchino fc Bvkovl l|200lt) . while Src 11 and 12 correspond to 1SAX 



parameter is < 20" in 97% of the cases, and therefore 
indicates reliably if a source is point-like or not. 

2.2.2. Sources identification 

In order to investigate the nature of the hard X-ray sources 
in IC443, we have carried out an extensive set of cross 
correlations with existing catalogues and digital data. In 
particular, we have used all the catalogues present in 
the W3Browse facility of Goddard Space Flight Center 
of NASA, and in the SIMBAD facility of the Centre de 
Donnees astronomiques de Strasbourg. We have used a 



standard search radius of 2' for all the catalogues except 
astrometric optical catalogues such as GSC, USNO and 
Tycho, were we have used 10". The list of identification is 
reported in Table El 

To verify if the hard X-ray sources are due to the in- 
teraction between IC443 and the molecular cloud w hich 
lies in foregr ound llBurton et alJJl98sl iR.ho et al.l200l|) . as 
suggested bv lBocchino fc Bvkovl l)2000|) . we have exploited 
the high spatial resolution of the PN data for a compar- 
ison with high-resolution optical and infrared data. We 
have used the Digital Sky Survey and the near-in frared 
data of the Two Micron All Sky Survey (2MASS, ICiitril 
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N RA Dec WCORR a FWHM cts s" 16 T cxp Note 

J2000 arcsec xlO" 3 sec 



1 


6 


15 30.1 


22 


49 


16.7 


5.8 


12.0 ± 7.2 


36.7 


± 


3.6 


7107 


s 


2 


6 


15 32.1 


22 


39 


33.0 


0.7 


11.9 ± 2.8 


6.3 


± 


1.3 


9335 




3 


6 


15 54.6 


22 


50 


55.6 


1.3 


11.3 ± 1.8 


5.7 


± 


1.2 


9994 




4 


6 


16 27.1 


22 


17 


18.2 


1.5 


10.5 ± 1.4 


2.8 


± 


0.6 


20441 




5 


6 


17 05.5 


22 


21 


29.7 


11.0 


20.5 ± 1.2 


65.8 


± 


4.4 


16301 


S, Plerion 


6 


6 


17 06.3 


22 


26 


51.1 


0.5 


12.5 ± 3.8 


4.8 


± 


0.9 


18356 




7 


6 


17 29.6 


22 


44 


53.9 


0.7 


11.0 ± 2.6 


7.0 


± 


1.5 


9398 


S 


8 


6 


17 33.9 


22 


25 


23.4 


2.1 


11.7 ± 1.3 


5.8 


± 


0.9 


22860 




9 


6 


17 40.2 


22 


52 


10.2 


0.2 


53.0 ± 22.0 


13.1 


± 


2.3 


4430 




10 


6 


18 01.5 


22 


22 


33.7 


1.3 


11.3 ± 1.7 


7.2 


± 


1.2 


13382 


s 


11 


6 


18 04.3 


22 


27 


32.8 


2.9 


20.6 ± 2.2 


24.7 


± 


3.2 


12463 


s 


12 


6 


18 06.4 


22 


28 


32.8 


3.3 


12.2 ± 1.0 


23.4 


± 


2.3 


11471 


s 



a The maximum correlation coefficient is a measure of source significance. 

b The count rate is given in the 2-10 keV band for comparis on with other publ i shed r esults. 
The flux was corrected for the PSF following the suggestion of lValtchanov et etl.l (l200lT) . 10 -3 
cnt s _1 is corresponding to 7.8 x 10~ 15 erg cm -2 s _1 adopting a power-law emission model 
with 7 = 2 and interstellar absorption Nh= 7 x 10 21 cm -2 . 

c "S" marks sources with a counterpart in the soft band. 

Table 3. Identification of PN detected sources with catalogues. 



N 


Source name 


Name 


Catalog 


Offset 
arcsec 




Remarks 


1 


XMMU J061530.1+224917 


U061530.75+224910.6 


USNO 


9 




B=18.8,R=15.8 


2 


XMMU J061532. 1+223933 












3 


XMMU J061554.6+225055 


06128+2250 


IRASPSC 


84 


/l2,25,6l 


,,100 = 3.1, 0.9, < 1.1, < 8.2 Jy 


4 


XMMU J061627. 1+221718 












5 


XMMU J061705.5+222129 


U061705.98+222120.0 
1WGA J0617.1+2221 
H061706.08+222131.3 


USNO 
WGACAT 
ROSHRI 


11 
6 
6 




B=16.7,R=15.9 






1SAX J0617.1+2221 






Src A 


in Bocchino & Bvkov (2000) 


6 


XMMU J061706.3+222651 










7 


XMMU J061729.6+224453 


U061729.57+224452.8 


USNO 


1 




B=15.9,R=15.6 


8 


XMMU J061733.9+222523 


TXS 0614+224 


TEXAS' 1 


34 




PosErr < 1" 


9 


XMMU J061740.2+225210 












10 


XMMU J061801.5+222233 


H061801. 60+222229.8 
RICH88 


ROSHRI 
DIXON 6 


4 
74 






11 


XMMU J061804.3+222732 


1SAX J0618.0+2227 c 






Src B 


in Bocchino & Bvkov C2000~) 


12 


XMMU J061806.4+222832 


1SAX J0618.0+2227 c 






Src B 


in Bocchino & Bvkov C2000~) 



° The Texas Survey of Radio Sources at 365 MHz oflDouglas et alJ (Il996ft . 
6 The Master List of Radio Sources (Version 43) of lDixorJ dl970l) . 

c The BeppoSAX source 1SAX J0618.0+2227 is resolved in two sources by our XMM-Newton observation. 



Il998l) . which has a resolution of ~ 3". In particular, we 
have considered the K s (2.32 /mi) band of 2MASS, for 
which most of the emission can be explaine d by molec- 
ular H 2 rovibrational lines l|Rho et alJl200l|) . a classical 
indicator of a shocked molecular cloud. Fig. shows the 
correspondence between the infrared, optical and X-ray 
emission of the hard X-ray sources located in the region 
of interaction between the SNR and the molecular cloud. 



Table 4. Definition of hardness ratios. 



Name 


Si 


Hi 


HRi 


0.2-0.5 


0.5-2.0 


HR 2 


0.5-2.0 


2.0-4.5 


HR 3 


2.0-4.5 


4.5-10.0 



2.2.3. Hardness ratios of detected sources 
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Fig. 3. Infrared/optical/radio images of the putative region of interaction with a molecular cloud. Left: optical DSS 
image. Center: 2MASS K s (2.32 /mi) band image. Right VLA NVSS radio survey at 1.4 GHz. All images are 18' size 
and we have superimposed the position of hard X-ray sources using a circle of 20" radius. 



In order to characterize the spectral properties of the de- 
tected source, we have computed three different and sta- 
tistically independent hardness ratios (HRs) of the form 

HR i = §^; i=l,2,3 (1) 

where the spectral bands are d efined in Tab ic IH The 
definition is the same as adopted by Ha singer et a l. ( 2001) 
to classify the sources detected in a deep XMM-Newton 
observation of the Lockman hole. The hardness ratios of 
the detected sources are displayed in Fig. in which we 
have also reported, for comparison, the HRs of a sample 
of regions of the thermal shell (Thl - Th4), and the four 
region s of the plerion nebula defined in lBocchino fc Bvkovl 
( 200l|. The latter can be described by a pure non-thermal 
power law spectrum with a slope steepening toward the 
outer regions (i.e. going from the core t o the "pl3" region, 
see Table 2 of lBocchino fc BvkovlEoOll) . 

Most the sources, except source number 7, clearly have 
a non-thermal origin, and are located in a well-defined re- 
gion of the diagram, near "pl3" , i.e. the most external ple- 
rion region which has a power-law spectrum with a photon 
index of 2.3. Source 1, 11 and 12 are located in a different 
group, near the regions "psr" and "pll" , thus implying an 
intrinsically harder and/or more absorbed spectrum. 

2.2.4. Spectral analysis of the sources 

We have chosen to extract the spectra of the bright- 
est sources which correlate with the infrared filaments, 
namely Src 10, Src 11 and 12 (Fig.|2J). A background re- 
gion outside the IC443 bright thermal emission has been 
identified. We have also identified the Th4 region at ~ 1' 
from Src 12 to be used as background region for the anal- 
ysis of Src 11 and 12 which are surrounded by substantial 
thermal emission from the SNR. In all the cases the back- 
ground regions have been corrected for vignetting and area 
effects prior to subtraction from source spectra, using the 



Table 5. Summary of PN spectral fitting results. 



Name 


N H 
xlO 21 


7 


flux" 
erg cm -2 


s- 1 


x'/dof 


SrclO 
Src 11 
Src 12 


1.6t^ 

8-li?;* 


l-86t^ 

i-48l° ; 2 ° 

9 ni +0.27 
^• O1 -0.24 


5.0 x 10" 
4.0 x 10" 
2.3 x 10" 


-14 
-13 
-13 


14/17 

97/88 
64/67 



" Absorbed power-law flux in the 0.5-10 keV band. 



ancillary response files (arf) computed with the SAS task 
ARFGEN. As a final step, we have summed the spectra of 
MOS1 and MOS2, averaging the arf files, and we have 
used the response matrix of MOS1 for the analysis of the 
summed spectrum. The fittings have been performed on 
the PN and MOS-summed spectrum for each source. 

Src 11 and 12 were co nfused into a single sour ce in 
a BeppoSAX observation ijBocchino fc Bvkovl |2000(L but 
are clearly separated by XMM-Newton into one extended 
source (n. 11) plus a point-like source (n. 12, Fig- EJ ■ Given 
the vicinity of these two sources to the diffuse thermal 
emission, we have analyzed their spectra using the Th4 
region as background. We have extracted the spectra from 
a circle with radius of 16" and 24" for Src 12 and 11, 
respectively. We have applied the PSF correction to the 
effective area used in the fits of Src 12 spectrum only. The 
spectra of these two sources are shown in Fig. Altogether 
with their best-fit model, while the best-fit parameters 
and their uncertainties are reported in Table A fit to a 
combination of thermal models (not reported) shows that 
the thermal nature of the spectra of these sources can be 
ruled out safely. 

The spectrum of Src 10 is also non-thermal. In spite 
of the higher uncertainties associated with its best-fit pa- 
rameters (this source is sensibly fainter than the Src 11 
and 12, see Table |SJ, we can say this source is character- 
ized by a spectrum of intermediate slope between Src 11 
and 12, and a significantly lower absorption. 
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Src1 1 



1.00 




Powerlaw 



-0.6 -0.4 -0.2 
HR2 (0.5-2 vs 2-4.5) 



0.0 




-0.8 -0.6 -0.4 -0.2 0.0 0.2 
HR3 (2-4.5 vs 4.5-1 keV) 

Fig. 4. Hardness ratios of the detected sources (1-12). We 
have also plotted for comparison the HRs of a sample of 
thermal (some knot of the shell diffuse emission, thl — thA) 
and non-thermal (the plerion nebula, pll — psr) sources. 
The expected HR values for a power-law model are also 
shown as solid lines (constant N H = 10 20 , 10 21 , 10 22 , 10 23 
cm~ 2 , from bottom to top) and dotted lines (constant 
7 = 1,2,3 from right to left), 
bars are shown. 



Some representative error 



3. Discussion 

The XMM-Newton observations of IC443 presented above 
have unprecedented sensitivity in 2-10 keV regime. The 
nature of the detected localized sources within the bound- 
ary of IC443 SNR should be a subject of a careful inspec- 
tion because some of them could be extragalactic sources 
unrelated to the SNR. Using the sensitivity map of the de- 
tection run,we have computed the log(N)-log(S) relation 
in 2-10 keV band of our IC443 observations, and we have 
compared it with the same function computed for the set 
of 10 XMM Ga l actic P lane Survey (GPS) observations of 
iBocchino et ail (l2002h . The GPS log(N)-log(S) is system- 
atically lower than the relation we have found in IC443 by 
a factor of 2 in the range 10 -3 — 10~ 2 cnt s _1 , and even 




Srd 2 



Powerlaw 




Fig. 5. PN+MOS spectra of Src 11 and 12 with power- 
law best-fit model. The best-fit parameters are reported 
in Table El 



more for brighter sources, thus indicating a clear excess of 
bright sources in the IC443 field of view. It seems there- 
fore likely that some of the detected sources are expected 
to be really related to the SNR on statistical grounds. We 
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Fig. 6. Confidence level contours (68%, 90% and 99%) 
in the Nh— 7 parameter space for Src 10, 11 and 12. We 
have also plotted the expected range of Nh for IC443 and 
G189.6 +3.3 overlapping region, for IC443 only (both ac- 
cording ^ao^^^^^chenbac^^^), and the total galac- 
tic Nh in the direction of IC443. 



have therefore investigated the possible nature of compact 
X-ray sources inside IC443. 

The hardness ratio distributions (Fig. allow us to 
conclude that only Src 7 is most likely of galactic nature, 
probably an X-ray counterpart of an active star. As for the 
other sources, the HR diagram does not make it possible 
to tell if they are galactic or extragalactic. The results of 
the spectral analysis on some of the sources give us ad- 
ditional informations. In Fig. |HJ we report the contour of 
confidence levels for the interstellar absorption and the 
power-law spectral index of Src 10, 11 and 12, along with 
some expected values for Nh- We can conclude that Src 
10 and 11 have the highest probability to be of galactic 
origin, while for Src 12 there is a large overlap between 
the allowed range of Nh and values well above the total 
galactic Nh toward this direc tion. However, we note that 
lAsaoka fe Aschenbachl l|1994|) report a total absorption of 
~ 10 22 cm" 2 for regions inside the shells of IC443 and the 
overlapping G189.6+3.3 SNR, as the region we are consid- 
ering here. We also note that, in addition to the absorption 
along the line of sight, there could be a considerable con- 
tribution to Nh fro m the molecular cloud. For instance, 
iRichter et alJ l|l995h have measured an extinction in the 
2.2 /xm band of ~ 1.5 in this region, which corresponds to 
a value of Nh^ 3 x 10 22 cm" 2 . Therefore, after all, there 
is no strong evidence even for Src 12 to be extragalactic. 

Having established that there may be a relation be- 
tween the sources and IC443, we now discuss a differ- 
ent scenario which may explain the association. There 
are three basic classes of compact hard X-ray sources re- 
lated to the SNR. One of these, the existence of a pulsar 
wind nebula inside IC443, has already been discussed by 



lOlbert et alJ ll200ll) and lBocchino fc Bvkovl ll200ll) for Src 

5, and therefore we will not discuss it here. 

3.1. Shocked molecular clump 

Shocked molecular clumps could represent another class 
of compact hard sources. As was discussed earlier, IC443 
is known to be interacting with a molecular cloud and 
the southern part of the remnant reveals a clumped and 
knotty structure of the IR molecular e mission down to 
arcsecond s cales (e.g.lRichter et al.lll995|) . Recent ISO ob- 
servations l|Cesarskv et alJll999|) of the southern region 
are consistent with that expected from shocked molecular 
clump models. This is also shown by Fig. which shows 
the optical, infrared, and radio images of the region in 
which the interaction with the cloud is occurring. In par- 
ticular, the absence of bright optical filaments is due to the 
ve ry low shock velocity and very high density confirmed 
bv lRichter et alJ l|l995T) . while the strong emission in the 
2.2 /im image is primarily due to the 1-0 S'(l) line of H 2 at 
2.122 /im, another strong indicator of slow shock in molec- 
ular clouds. It is remarkable that 5 out of 11 (excluding 
the plerion) compact X-ray sources are located close to 
this infrared-filament network. 

Modeling of nonthermal emission from an SNR inter- 
ac ting with a molecula r cloud has been recently performed 
bv lBvkov et a l. ( 2000). It has been shown that hard X-ray 
and 7-ray emission structure should consist of an extended 
shell-like feature related to a radiative shock, and local- 
ized sources corres ponding to shocked molecular clumps. 
iBvkov et, alJ l|2000l) predicted a hard spectrum with a pho- 
ton index below 2.0. Because of the hard spectrum and 
heavy absorption, localized spots of a few arcmin size 
would be seen only in hard X-rays. The emission orig- 
inates from shock-accelerated electrons of energy below 

I GeV. The predicted flux density is ~ 6 x 10" 5 keV 
cm" 2 s _1 keV" 1 at 10 keV for a 30 km s" 1 shock in 
a clump of half parsec radius and density of 10 4 cm" 3 . 
The predicted count rates are c onsistent with the observe d 
rate of 1SAX J0618.0+2227 iBocchino fc Bvkovl Eoooh . 
Recent ASCA study of 7 Cygni, a shell-type SNR inter- 
acting with a molecular cloud, revealed several clumps of 
hard X -ray emission dominati ng the SNR emission above 
4 keV llUchivama et al.ll2002j) . While somewhat brighter 
(~ 4.5 x 10 -12 erg cm" 2 s" 1 ) and more extended, the 
clumps resemble the spectra and surface brightness of the 
source 1SAX J0618. 0+2227 inside IC 443. 

Therefore, 1SAX J0618.0+2227, located near bright 
spots of molecular hydrogen emission, could be a shocked 
molecular clump. However, with the XMM-Newton PN 
camera the source 1SAX J0618. 0+2227 is resolved into 
two different sources, Src 11 and 12 (Fig. |2J). Src 12 ap- 
parently has a soft spectrum and a size of 12" FWHM, 
which makes the clump interpretation unlikely. The Src 

II has a hard spectrum, but still only a size of 20" 
FWHM (~ 4 x 10 17 cm at 1.5 kpc). The shocked molecular 
clump interpretation of Src 11 would require a substan- 
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tially higher shock velocity in the clump than 30 km s _1 
to match the observed flux. The energetic requirements 
are rather stringent for such a small scale shock. 

3.2. Fast-moving isolated ejecta fragments 

The third class of possible compact hard X-ray sources 
related to the SNR are the isolated ejecta fragments. 
Multiwavelength studies of the structure of supernova 
remnants have revealed the presence of fine-scale struc- 
ture with rich emission spectra. There are metal-rich frag- 
ments of SN ejecta interacting with surrounding media 
as well as circumstellar mass-loss matter. Some X-ray 
knots, e.g. th e silicon-rich shrapnel, w ere o bserved in the 
Vela S NR by lAschenbach et alJ l|l995l) and iMivata et ail 
ll200ll) and also in t he Ty cho SNR with XMM-Newton by 
iDecourchelle et ail l)200ll) . Supersonic motion of the knots 
in the ambient medium will result in a bow shock/knot- 
shock structure creation. It was shown that an individual 
knot propagating through a molecular cloud would be ob- 
servable with XMM-N ewton and Chandra from a few kpc 
distance l)Bvkovl2002|) . Such a know would have a hard X- 
ray spectrum of photon index < 1.5 in 1-10 keV regime, 
with prominent lines of some metals e.g. Si, S, Ar, Ca, Fe 
etc. 

In order to verify the knot hypothesis for Src 11 and 
Src 12, we investigated the presence of emission lines in 
the X-ray spectrum. We added a gaussian line of variable 
width in the range 0-200 eV to the power-law model and 
we fitted the spectra, thus deriving the \ 2 contour levels 
for the line centroid and flux shown in Fig. We did 
not detect any emission lines at the 99% confidence level, 
even though there is excess flux at ~ 1.8 keV, 3.4-3.6 keV 
and 5.1 keV only for Src 11. Fig. reports the upper- 
limit to the line flux at different centroid energies. We 
have verified that the excess flux cannot be explained by 
known background emission lines. 

We applied the supersonic knot model of lBvkovlll2002l) 
to Src 1 1 . X-ray line emission was simulated from a metal 
rich ejecta fragment of 3 x 10 17 cm radius (~ 15" at 1.5 
kpc) and a mass of ~ 10~ 2 M Q interacting with a molec- 
ular cloud. The simulated fragment mass was dominated 
by oxygen with a minor mass contribution (~ 10 _4 Mq) of 
other SN nucleosynthesis products like Si, S, Ar, Ca, Fe. 
The knot travels through an inter-clump medium of num- 
ber density < 100 cm" 3 at a velocity of ~ 500 km s _1 . 
Under these conditions, the model yields a 1.8 keV Si line 
flux of > 10 -6 ph cm -2 s _1 (if the electron diffusion co- 
efficient is < 10 18 cm 2 s _1 ), to be compared with our 
upper-limit of 4 x 10 -6 ph cm -2 s _1 . The optical depth of 
the 1.8 keV Si line (Si VI and higher ionization states) is 
~ 6 in that case, implying a saturation effect, i.e. increas- 
ing the Si mass would not increase the observed flux, if the 
knot radius is the same. The depth is somewhat lower for 
heavier elements where the saturation becomes improtant 
only above 10 _4 Mq. A higher knot velocity (or lower dif- 
fusion coefficients) would produce higher Si line fluxes, as 
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Fig. 7. Confidence level contours (68%, 90% and 99%) in 
the line centroid vs. line flux parameter space for Src 11 
(solid lines) and 12 (dotted lines). We have used a model 
with a single gaussian line of variable width plus a power- 
law continuum. There is no firm detection of any line; the 
99% contour levels (the uppermost for each source) define 
our estimate of the upper-limit to the line flux. 

so would a bigger knot. However, the apparent size of Src 
11 would point against knots having sizes sensibly differ- 
ent from those used in our simulation. The knot position, 
given the SNR age above 10 kyr, is constraining the knot 
velocity. Thus the electron diffusion coefficient is the most 
important parameter, which may be varied to match more 
accurate measurement of the 1.8 Si emission line flux. 

We note that the lifetime of the ejecta fragment in 
the inter-clump medium is rather short, being below 100 
yr. This implies that the fragment should spent most of 
its time in a tenuous gas of density < 0.1 cm~ 3 of the 
SNR interior to survive during the lifetim e of IC443, whic h 
is expected to be about 30,000 yr (e.g. IChevalierlll999l) . 
Indeed, the Src 11 is located just on the border of the 
molecular cloud, which suggests that the interaction with 
denser material has just started. Note that the [Ar III] 
bright knots of velocities ~ 500 km s" 1 (but with a some- 
w hat smaller size) observed in the Crab nebula with HST 
bv lSchaller fe Fesenl l|2002j) could represent a similar phe- 
nomenon 

The definite test to distinguish between the two in- 
terpretations, namely a fast moving metal-rich knot or a 
shocked clump, would be an X-ray spectral lines study 
because the shocked molecular clump emission would be 
dominated by a hard continuum in contrast to the metal- 
rich knot case. High-quality line spectra of the Src 11 
would allow the X-ray study of SN nucleosynthesis. 

4. Summary and conclusions 

We have detected 12 compact hard X-ray sources in the 
field of the large SNR IC443, more than the expected num- 
ber of galactic and extragalactic sources in this field, thus 
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suggesting an association with the SNR for some of them. 
We have carried out a spectral analysis of the sources, 
by fitting the X-ray spectra of the brightest ones, and 
by computing X-ray hardness ratios. We have found that 
all sources but one have a clearly non-thermal spectrum. 
We have noted that 6 (including the plcrion) out of 12 
sources are concentrated in a small (~ 15') area of the 
SNR, which was previously known to be a region of inter- 
action between the SNR shock and a molecular cloud. We 
suggest that some of the sources may indeed be related 
to the SNR itself, and we review the possible formation 
mechanisms of formation of isolated hard X-ray sources in 
SNR shell. We have considered both electron acceleration 
in molecular clump shocks, and the emission from ejecta 
fragment as an alternative model. The latter is suggested 
by a modest evidence of a Si line in the spectrum of one 
of the sources. A deeper X-ray spectroscopy observation 
of the sources in IC443 and more X-ray observations of 
other SNRs interacting with molecular clouds are urgently 
needed to shed more light on this intriguing new class of 
hard X-ray sources. 
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